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We derive a general microscopic expression for the non-linear diamagnetic current in a clean 
superconductor-insulator-normal metal structure with an arbitrary interface transmission. In the 
absence of electron-electron interactions in the normal metal the diamagnetic response increases 
monotonously with decreasing temperature showing no sign of paramagnetic reentrance down to 
T = 0. We also analyze the magnetic breakdown of proximity-induced Meissner screening. We 
demonstrate that the magnetic breakdown field should be strongly suppressed in the limit of small 
interface transmissions while the linear diamagnetic current does not depend on the transmission of 
the insulating barrier at low enough temperatures. 



I. INTRODUCTION 

Mesoscopic hybrid structures composed of supercon- 
ducting (S) and normal (N) metals demonstrate a reach 
variety of intriguing physical phenomena. Many of them 
have recently been extensively investigated - both ex- 
perimentally and theoretically - and received adequate 
interpretation within the framework of the quasiclassi- 
cal theory of superconductivity, see, e.g., Refs. 0J| for a 
recent review. 

One of the remaining puzzles in the field concerns the 
paramagnetic reentrance phenomenon observed in silver- 
niobiumd'Q and gold-niobiumQ proximity cylinders in the 
limit of very low temperatures. This reentrance behav- 
ior is in aj-|Clear contradiction with earlier theoretical 
predictipnsS as well as with the results of more recent 
studiesQ demonstrating that diamagnetism in SN prox- 
imity cylinders should progress monotonously with de- 
creasing temperature. Several theoretical explanations-of 
this reentrance phenomenon have been put forwardtjLltLll 
However, a detailed comparison with experiments carried 
out in Refs. ^,[| indicates that none of these explanations 
is able to correctly reproduce the absolute value or the 
temperature dependence of the observed paramagnetic 
reentrance effect. One of thej-epcplanations has also been 
subject to theoretical debateEil. 

It is important to emphasize that, while the authorsiB 
invoke additional assumptions about the form of electron- 
electron interactions in the N-mctal, the worku oper- 
ates with the standard model for SN proximity systems 
which assumes no interaction between electrons in the 
normal layer. This model is usually very well described 
by means of the standard approach based on the qua- 
siclassical .Eilenberger equations. In order to test the 
conjecturecl, that paramagnetic reentrance could be miss- 
ing within the standard quasiclassical formalism because 
of an additional contribution from the low energy "glanc- 
ing" states, one needs to go beyond quasiclassics and ana- 
lyze the problem within a more general approach. Below 
we will use the microscopic Gor'kov equations and eval- 



uate the screening current by constructing the full Green 
functions for the problem in question. We will then com- 
pare our result with that derived from the Eilenberger 
equations. This comparison confirms the applicability of 
the standard quasiclassical technique and demonstrates 
that no paramagnetic reentrance can occur in SN prox- 
imity cylinders in the absence of electron-electron inter- 
actions in the N-layer. 

Another interesting property of the same structures is 
the so-called magnetic breakdown of proximity-induced 
Meissner screening. This effect occurs at a certain value 
of the external magnetic field H^, above which the dia- 
magnetic expulsion of the field from the norma] layer 
turns out to be energetically unfavorableoE3E3 although 
the proximity effect itself is not yet suppressed. In the 
absence of impurities in the N-layer and for the case of a 
perfectly transparent interface between S- and N-metals 
the breakdown field H^(T) was evaluated in Ref. 
Since at sufficiently low temperatures the linear current 
response of SN proximity cylinders [does not depend on 
the transmission of the SNjaterfacaiJ one could assume 
that at T — ► the resultta for should also apply 
at arbitrary transmissions. However, this turns out not 
to be the case. As we will show below, the breakdown 
field scales linearly with the interface transmission and, 
hence, should be strongly suppressed in the limit of 
small transmissions. This effect makes it possible to ex- 
tract information about the interface transmission from 
the measurements of the magnetic breakdown field H^. 

Our paper is organized as follows. In section II we will 
present a general derivation of the non-linear screening 
current in superconductor-normal metal proximity sys- 
tems with arbitrary interface transmissions. This deriva- 
tion is carried out on the basis of the microscopic Gor'kov 
equations and does not involve energy-integrated quasi- 
classical Eilenberger functions. The effect of magnetic 
breakdown of proximity-induced screening is analyzed in 
section III. Applicability of the quasiclassical approach 
for the problem in question is discussed in section IV. 
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II. NONLINEAR DIAMAGNETIC RESPONSE 

We will consider a clean SIN structure and assume it to 
be uniform along the directions parallel to the interfaces 
(coordinates y and z). The normal metallic layer (N) 
and the bulk superconductor (S) are located respectively 
at < x < d and x < 0. Our analysis will be based 
on the Gor'kov eouations of the microscopic theory of 
superconductivity£3. After the Fourier transformation 
of the normal (G) and anomalous (F + ) Green functions 
with respect to y and z, 



the Gor'kov equations take the following standard form 



iuj n — H 
A*(x) 



iu) n + H c 



G Un (x,x' ,k\\] 
Fw n {x,x',k n ) 



S(x — x' 




(1) 

Here u> n = (2n + l)irT is the Matsubara frequency, and 
A(x) is the superconducting order parameter. Below we 
will choose A (a;) = A for x < and A(x) =0 otherwise. 
The Hamiltonian H in Eq. (nT) is defined as 



H 



1 d 2 
2m dx 2 



k 2 

2m 



e F + V(x), 



(2) 



where fey = fey — -Aii(x), ep is Fermi energy and the 
term V(x) accounts for the boundary potentials. The 
Hamiltonian H c is obtained from H (|^) by inverting the 
sign of the electron charge e. 

In what follows we shall neglect the square term oc A 2 ^ 
in Eq. (^). This approximation applies in the range of 
magnetic fields H -C &okp/d, where kp is the Fermi 
momentum and <I>o = nc/e is the superconducting flux 
quantum. We will assume that the London penetration 
depth of the bulk superconductor is small and neglect 
the magnetic field inside the superconductor. The latter 
assumption allows to set A(x < 0) = 0. 

The next standard step is to decompose the Green 
functions into the product of quickly oscillating terms 
exp(±ik x x) and the two component envelope functions 
Tp±(x) changing at scales much longer than the Fermi 
wavelength. One has 



iuj n — H A(x) 
A* (x) iw n + H, 

a ±zkxx ( iu n - H± A(x) 



ip±(x)e 



where we defined k x = 4 / k 2 F — k 2 and 



H£ = Tiv x d x •Ai|(x)t>|| 



(3) 



(4) 



A convenient choice of the gauge for the problem in ques- 
tion is A = (Q, A(x),0). Below we will make use of 



this gauge aa.d proceed along the lines with our previ- 
ous analysistJ. 

Let us fix the coordinate x' within the normal metal. 
In this case for x < the general solution of Eq. (|l|) can 
be written in the form 



Gui n {x, x') 
F+(x,x') 



1 

—i 



e Ax/v *e lk * x si{x') 



(5) 



s 2 (x'). 



Here we restricted ourselves to small Matsubara frequen- 
cies ui n <C A. This is sufficient in the most interesting 
physical limit r< A and d ^> £o ~ vp/A which we only 
consider below. 

The general solution of Eq. (ffl) for x > reads 



G Un (x,tf) \_f G Un (x,x') 
F+Sx,x') ) \ 

~'h{x') 
'h(x') 



(6) 



(p^(x)e lk * x f 3 (x') + <p+{x)e 



— ik x 

— ik x 



Here the functions 

ip + (x 



( w n f X e ^v a i si 

cxp I x + i I —-^-Ayxijaxi 

\ v x Jo c kx 

1 

f-( x ) 



(p+(x) 



(7) 



obey the equations Hlw±(x) = and the first term in 
the right hand side of (||) with 

r /_ „n = _i_ [ ip+(x)ip-(x')e ik ^ x - x "> iix>x' 
Un{ ' ' v x \ ^{x) V+ {x')e lk ^ x '-^ if i < x' 

_ (8) 

represents the particular solution of Eq. for a; > 0. 
What remains is to determine six functions Si^(x') and 
fi-4(x') in Eqs. ^^)- This is done with the aid of the 
boundary conditions imposed on both sides of the normal 
layer. 

First we consider the SN interface. Matching the wave 
functions on both sides of this interface, respectively 
A\ exp(iki x x) + B\ exp(— ik\ x x) and A2 exp^A^x) + 
B2 exp(— ik2 X x), is performed in a standard manner (see, 
e.g., Ref. ph: 



A 2 = olA x + 0B X , B 2 = P*Ai + a*B x , 



a 



1, 



(9) 



where the reflection and transmission coefficients are 
equal to 



R 



D = 1 - R = 



1 



(10) 



Applying Eqs. (||) to the two-component vectors (G, F + ) 
we get four linear equations for six unknown functions in 
Eqs. (§,§). The remaining two equations are derived 
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at the interface between vacuum and the normal metal 
(x = d) where we assume complete specular reflection of 
the electrons. This boundary condition trivially yields 



G{d,x') = F+(d,x') = 0. 



(11) 



The above six linear equations uniquely define the Green 
function of our system. Technically the calculation is 
similar to that presented in Ref. 15, therefore we will 
not go into further details here. The resulting expression 
for the Green function in the N-metal (0 < x < d) takes 
the form 



Guj n {x, x, fcy ) — 



1 



zsinhx + fxjf sin (7) 



coshx 



$gcos( 7 ) 



X 



2u> n d 



2i 



^Y^-A(x)dx, 7 = 2k x d + ar 5 — 
ck T a 



(12) 



Of interest for us here is the current density in the y- 
direction. It reads 



D = 0) Eqs. (|4 
derived in Ref. 



( fl7|) reproduce the non-linear response 
Another important limit is that of a 
small external field. In this case one can linearize the 
function G © in d> and find 



ekpTip y--^ 



uj n >0 Jo 



n/2 



dd sin 3 9t 2 (9) coshxi 
cos 9 [t 2 {9) +sinh 2 xi] 



3/2 - 



(18) 

This result was obtained in Ref. [lj with the aid of a 
different approach. It is interesting that at low temper- 
atures T <C tvp/d the transmission and reflection coef- 
ficients drop out and for any nonzero t the result ( |l8| ) 
reduces to a simple formula 



3 = 



ekpVpff) 
6n 2 d ' 



(19) 



which was initially derivedQ for t = 1 . At higher temper- 
atures t does not anymore drop out of the final result. 
For instance, in the limit T 3> vf / 'd we obtain 



m = -t e 

171 u; n >0 J \h\<kF 



dk \\ 

— — fcyRe Guj n (x, A; || ). 



(13) 

Note, that the terms sin 7 and cos 7 in Eq. (O) are 
quickly oscillating functions of fey. Integp±ion over the 
momentum directions in (|l^) is equivalent!!^ to averaging 
over the angle 7. Performing this averaging we arrive at 
a spatially homogeneous screening current in the N-layer 

j = f— 2_] / d9 dip sin 2 cos i^Im G(%)- 

_ (14) 
The function G in Eq. (|lj) depends on the complex 
variable \ — Xi — *X2, where 



Xi 



and 



2uj n d 
vf cos^' 



\2 = 4> tan 9 cos ip (15) 



A{x)dx. (16) 



The function G{\) is 7r-periodic in \2- Within one period 
it is defined on a strip — 7r/2 < \2 < with the cut 
going from (0, — arcsini) to (0, arcsini). For such values 
of X2 one has 



G 



sinh x 



t 2 (9) +sinh 2 x 



(17) 



where t{9) = D(6)/(l + R(9)). Eqs. @-@ m % de- 
termine the non-linear diamagnetic response of a clean 
SIN system to an externally applied magnetic field H . 

It is easy to check that our result reduces to the already 
known ones in the corresponding limits. For instance, in 
the case of a perfectly transmitting SN interface (R = 1 — 



J 



ek F v F t cj) / t^\ ^ Td/VF 



2n 3 d \TdJ 



(20) 



where to is the value of t{9) at 9 = 0. In the case of low 
transmissions to -C 1 there exists an additional temper- 
ature interval toVF/d -C T <C vp/d, where the current 
shows a power-law dependence on temperature: 



. _ 7C,(Z)ek 2 F v F (j) /v F V 
3 _ 647r 4 d VTd/ 



ir/2 



d6> cos 2 sin 3 W 2 (6»). 



(21) 

It is important to emphasize that the screening cur- 
rent in the N-metal is uniform in space and, hence, is 
a nonlocal function of the vector potential. Combining 
this result with the Maxwell equations one easily arrives 
at the self-consistency condition for the "phase" 



^ = Hd 2 + ^ md \ 
e 6c 



(22) 



In the linear response regime the current density j 
can be expressed in the form j = — c 2 4>/^iieX 2 N (T)d. 
As soon as the effective length Xn(T) becomes small, 
Aat(T) -C d, the non-locality of the screening current 
turns out to be important. For instance, at T — * 
we have j — —3Hc/8nd, the magnetic field penetrating 
into the normal layer becomes over-screenecB, B(x) — 
3Hx/(2d) — H/2, and the average magnetization of the 
N-layer reads inM = -3H/A. 



III. MAGNETIC BREAKDOWN 

The above Meissner state of our hybrid structure 
is thermodynamically favorable only in relatively small 
magnetic fields. At higher values of H the diamagnetic 
screening current gets suppressed and the magnetic field 
can freely penetrate into the N-layer. This magnetic 
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breakdown effect was experimentally studied by Mjala 
and coworkcrsEj and was also addressed theoretically0'E3 
in the case of perfectly transparent SN interfaces. In this 
section we will consider the magnetic breakdown of the 
Meissner effect in clean proximity structures with an ar- 
bitrary transmission of the SN interface. We will closely 
follow the strategy adopted in Ref. [lj| 

The free energy F of the current state (normalized per 
unit surface) can be recovered by means of the standard 
procedure which amounts to integrating the current (|l4| ) 
over an effective "coupling constant" A 



F{4>) = -0/2e) / j{\<t>)d\ 



(23) 



In the presence of an external magnetic field we have to 
minimize the Gibbs free energy 



Q(T, H) = F{ct>) + ±- [ dx(d x A(x) - H) 
8tt Jo 



(24) 



In the limit of low external fields H < 7? m ; n the dia- 
magnetic solution described in the previous section is 
the only possible one. On the contrary, at high mag- 
netic fields H > i/max this solution cannot be real- 
ized. In this case another solution of the Maxwell equa- 
tions - which describes full penetration of the magnetic 
field into the N-layer - takes over. In the intermediate 
regime H m i n < H < -ff max the diamagnetic and non- 
magnetic solutions may coexist. The first order phase 
transition between these two states occurs in this inter- 
mediate regime at a certain value of the magnetic field 
H = iJb implying the magnetic breakdown of Meissner 
screening. Since this breakdown field can be suffi- 
ciently large, the full nonlinear expression for the dia- 
magnetic response, Eqs. (p^)-(|l7|), should be taken into 
account. 

First let us briefly consider the case of a perfectly trans- 
parent SN interface t = 1 and re-derive the resultsEj. 
The "supercooling" and "superheating" fields H m \ n and 
ff max can easily be estimated from Eq. (^||). Combin- 
ing this equation with the expression for the screening 
current j = —c 2 4>/8neX 2 ^(T)d one readily finds </> ~ 
3irHX 2 v (T)/§ - Since the diamagnetic solution is pos- 
sible only for small <\> -C 27r, one obtains 



$o/X 2 N (T). 



(25) 



For a non-magnetic solution the screening current is sup- 
pressed j w 0, in which case from Eq. ( |22] ) we find 
4> ~ irHd 2 /$o- The vanishing screening current implies 
strong dephasing of Andreev states in the N-layer. This 
dephasing effect occurs for cf> ^> 2tt. Hence, we get 



(26) 



Now let us evaluate the breakdown field H^. The 
Gibbs energy of a non-magnetic state Q nm is minimized 
by the equation d x A = H. Making use of Eqs. (Ill)- (ITT 



and (E3,E4) together with the inequality 4>^> 2tt one finds 



_ kpT x - 

yjim — / 

w„>0 



hi 



1 + exp - 



Au n d 

jJLVF 



lid/i, (27) 



where we denoted p = cos 9. 

The free energy of the diamagnetic state Qdm can also 
be easily established. Assuming Xn(T) Cfiwe obtain 



Q dm = 3H 2 d/32ir. 



(28) 



The breakdown field is derived from the condition 
Qdm = Qnm ■ In accordance with Ref. O one gets 



\/2$ 



-1nTdjv F 



nX N (0)d 

in the high temperature limit T ^> vpjd and 

$ 



6X N (0)d 



(29) 



(30) 



in the opposite limit of low temperatures T <C Vp/d. 
We observe that for Xn(T) -C d the breakdown field 
is parametrically larger (smaller) than H min (H max ), i.e. 
the phase transition indeed occurs in the intermediate 
regime iJ m ; n < H < H nmx . The condition Xn(T*) ~ 
d defines the temperature T* — (vp/2nd) ln(d/Ajv(0)) 
below which the above picture remains valid. At higher 
temperatures T > T* a continuous and-Deversible cross- 
over between the two states is cxpcctcdli3. 

Let us now turn to the case of a non-ideal SN interface 
t < 1. The estimate (^6j) for the field -ff m i n remains 
the same while (^5|) should be modified. Below we will 
estimate H max and obtain the expression for ifb in three 
different temperature intervals. 

High temperatures. In the limit r> vpjd Eqs. (|l4|)- 
© yield 



3 



zpvptfa 



2n 3 d \TdJ 



(VF\ iwTd/vp 



exp 



V F (, 



2nTd 



(31) 



Making use of this expression we evaluate the Gibbs en- 
ergy and find which is again defined by Eq. ( |29| ) 
multiplied by the factor to. Analogously, the estimate 
( p5| ) should now be multiplied by i 2 , and the whole pic- 
ture remains consistent at to > A at (T) j d. 

Intermediate temperatures. In the limit to -C 1 it is 
possible to realize the condition tovp/d <T« vp/d. In 
this intermediate temperature range we obtain 



t\6) COS X 2 



16ir 4 Td 2 Jo jo 
l + ip 



7T/2 

d9 I dip cos 2 9 sin 2 9 cos <p x 



1 — ip 



(32) 



where ip(z) — T'(z)/T(z) is the digamma function and 
vf cos 9 sin \i 



P : 



2nTd 



(33) 
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In the limit (f> <C TcL/vf we recover the result (|2lJ), while 
for <j) 3> Td/vp the screening current decreases with in- 
creasing </). For Td/vF C ^ 1 we get 



J 



ekpvp 
4iT 2 d(b 



\x t {jj)dfjb. 



(34) 



The estimate for H max is obtained by combining Eqs. 
( pl| ) and (p2|). In this regime -ff max turns out to be 
smaller than ( p5| ) by the factor ~ t^vp/Td. The free 
energy of a non-magnetic state is also established easily. 
Making use of the condition cf> ^S> 2ir and integrating the 
current (p2h over A in (p3|) we obtain 



F 



8n 2 d 



d^/i 2 t 2 (yLt) ln«, 



(35) 



where K = 7oM«f / (2TrTd) and In 70 ~ 0.577 is the Euler 
constant. The magnetic breakdown field is obtained 
by comparing Eqs. (p8|) and (|35|). One finds 



7rAAr(0)d 



<i/i/i 2 t 2 (/i) In k 



1/2 



(36) 



Thus, in the limit to <C 1 this field is strongly reduced as 
compared to the case of perfectly transparent SN inter- 
faces ©. 

Low temperatures. In the case T <C t$vp/d the dia- 
magnetic current takes the form 

ek 2 v F r /2 r /2 

j = I d9 I dip sin 2 9 cos 9 cos <pl(xi ■> *) • 

n d J Q J 

(37) 

Here I(x2, t(9)) is a 7r-periodic function of \2- For | X2 1 < 
7r/2 this function is defined as 

Vsin 2 X2 - t 2 



/ = X2 - signx2 6 (| sinx2| - arctan 

(38) 

where 0(x) is the Heaviside step function. The function 
I(X2,t) is depicted in Fig. 1 for different values of t. We 
observe that this function is linear in X2 oc only at 
cf> < t. For i < (4 < 1 we again recover the dependence 
(§)■ 

The estimates for if ma x and follow from the above 
expressions for the current. At T — > the field iJ max 
turns out to be by the factor ~ to smaller than ( gg ) 
and Hh is again defined by Eq. ( |36] ) with k ~ 1/to in 
this case. The above analysis remains consistent down 
to very small transmissions to ^ Ajv(0)/d. 

The dependence of on the transparency parameter 
to in the limit T — ► was also calculated numerically 
from Eqs. ( ^7| , ^8| ) and the result is presented in Fig. 
2. With a sufficient accuracy this dependence can be 
approximated by a simple formula 



H h (t )/H h (l) ~ to. 



(39) 




FIG. 1: The function Jfoa) for different values of t = 
0.2,0.5,0.8,0.99 (left to right). The initial slope remains the 
same for all curves. 




FIG. 2: The magnetic breakdown field as a function of to 
evaluated at T — for the angle-dependent transmission of 
the form @. 



The above results for H},(to) can be used to extract 
information about the effective transmission of the SN 
interface from the measurements of the breakdown field 
H\y. A detailed comparison between theoretical predic- 
tions (derived at t = 1) and experimental data for Hb was 
carried out in Refs. [13 18 It was found in both papers 
that for comparatively clean samples the temperature de- 
pendence of the breakdown field is quite well described 
by a theory derived in the clean limit. At the same time 
the amplitude of turned out to be somewhat smaller 
(typically by a factor ~ 0.3 -j- 0.6) than that predicted in 
the ballistic case. One can attribute this discrepancy to 
the effect of electron scattering on non-magnetic impuri- 
ties in the N-layerea It follows from the above analysis 
that another possible reason for this suppression of Hy, is 
a non-perfect transmission of the SN interface. 

For instance, one can argue that the difference between 
theory (30) and experiment by a factor « 0.56 found in 
Ref. 13 at low T could be due to the effect of the inter- 
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face transmission. Then from Fig. 2 one would recover 
t w 0.55. In the most relevant case of a thin effective 
potential barrier at the SN interface the angle-dependent 
transmission is defined as 



t{9) 



to cos 2 



(40) 



1 - t sin^ 

in which case the above estimate for to will translate into 
the value of the minimum interface reflectivity R m in — 
0.29. Hence, one can conclude that SN interfaces in 
the experimentsEj were indeed highly transmittive as it 
was already assumed before by a number of authors. Of 
course, other factors, such as electron scattering on im- 
purities in the N-layer, differences in geometry and cor- 
rections due to the term A 2 neglected in Eq. (Q), can 
cause additional minor discrepancies 2 !] between theoreti- 
cal and experimental values of Hh and can slightly mod- 
ify our estimate for i? m ; n . However, the above conclusion 
can hardly be affected by these minor modifications. In 
fact, the estimate for i? m i n could only increase if electron 
scattering on impurities is taken into account in addition 
to the effect discussed here. 



IV. APPLICABILITY OF QUASICLASSICS 
AND ABSENCE OF PARAMAGNETIC 
REENTRANCE 

Finally let us briefly address the applicability of the 
quasiclassical Eilenberger approach for the problem in 
question. Solving the Eilenberger equations in the S- 
and N-metals and matching these solutions at the SN in- 
terface with the aid of the Zaitsev boundary conditions 
we have evaluated the screening current in the N-layer of 
our proximity system. The calculation is to much extent 
analogous to that presented in Appendix A of Ref. [If], 
therefore we will not go into details here. Although the 
result of this calculation does not exactly match with our 
general expression ( |l3|), it turns out to be identical to one 
presented by Eqs. (|14|)-([l7|). The latter result was ob- 
tained from ( |l3| ) by means of averaging over all possible 
directions of the Fermi momentum or, equivalently, over 
the angle 7. Similarly to Ref. [l^ the difference between 
the expressions for the current before and after this aver- 
aging is small in the parameter 1/kpd. We also note that 
in the case of highly transmitting interfaces, R — > 0, no 
averaging over 7 is needed and Eq. ( |l3| ) coinddes exactly 
with the result of the quasiclassical analysisu. 

The above calculation was carried out in the limit of 
large spacial dimensions L y and L z of our system respec- 
tively in y and z directions. For finite L VtZ our results 
for the screening current acquire a small correction. Its 
magnitude can easily be estimated by rewriting the inte- 
gral over the momentum directions in Eq. ( |l3| ) as a sum 
over the conducting channels. Using further the Eulcr- 
Maclaurin summation formula 



1 



F(a) 



71=1 



F(a + n) 



F(x)dx - ^F'(a) (41) 



we observe that the correction to our results due to finite 
L yz is small as 1/a/AT, where Af ~ k 2 F L y L z is the effec- 
tive number of conducting channels in the x-direction. 

The parameters 1/kpd and 1/V77, therefore, control 
the accuracy of the quasiclassical Eilenberger approach 
in our problem. Both these parameters are very small 
since they contain the ratio between the Fermi wave- 
length and at least one of the system dimensions. For 
instancep- for the systems studied in Refs. p|,E|j5|Jl^ we 
estimated l/k F d < 10~ 4 and 1/VJ7 < 10 _e \ Hence, 
one can conclude that the Eilenberger formalism remains 
very accurate in this case and no significant physics is 
missing within this formalism at any temperature and at 
any interface transmission. 

The above estimates also indicate that within our 
model no paramagnetic reentrance effect can be expected 
in SN proximity cylinders. Indeed in the limit L y d 
the difference between the slab and cylinder geometries 
is negligible and our results will be directly applicable 
to the latter geometry if we identify L y with the cir- 
cumference of the cylinder. Then the correction from all 
electron states within the energy interval ~ A from the 
Fermi energy (obviously this correction also includes the 
contribution of the low energy glancing stateso) is small 
as - 

Although the conclusion about the absence of param- 
agnetic reentrance was obtained here only in the clean 
limit and for specularly reflecting interfaces, this conclu- 
sion should apply in the presence of impurities as well. 
In fact, a purely ballistic situation appears to be most 
favorable for the low energy glancing states. As soon as 
a finite electron mean free path I is introduced a proxim- 
ity induced miniaap ~ VpTXixnQ ,l/d ) will develop in 
the normal metalE3 and there will be no glancing states 
in the system at all. 

We also note that our present consideration does not 
include the contribution of the electron states with en- 
ergies smaller than ep — A. This contribution to the 
current (which has nothing to do with superconductivity 
and glancing states) can be roughly estimated with the 
aid of the standard results for the persistent currents in 
normal metallic cylinders, see, e.g., Ref. |23|. Making use 
of these results, in the limit of small magnetic fields one 

can assume Ipc ~ e(vp / ' L y )ky, 2 d 1 ^ 2 L z cf). The relative 
magnitude of this contribution is also much smaller than 
the proximity induced diamagnetic current studied here. 
At T — ► and for relatively high interface transmissions 
one finds Ipc/jdL, ~ [d/kpL 2 ] 1 / 2 . For typical experi- 
mental parametersflBi this ratio is of order ~ 10~ 4 , i.e. 
in this case no noticeable correction to our results could 
be expected either. 

In summary, with the aid of the microscopic 
Gor'kov equations we have derived a general expres- 
sion for the non-linear diamagnetic current in a clean 
superconductor-insulator-normal metal structure with an 
arbitrary interface transmission. We have demonstrated 
that at low temperatures the magnetic breakdown field 
H\j is suppressed linearly with the transmission of the in- 
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sulating barrier. This observation enables one to obtain 
information about the quality of inter-metallic interfaces 
from measurements of H\>. We have also compared our 
results with those derived within the quasiclassical Eilcn- 
berger approach and found the latter approach to be ex- 
tremely accurate for the systems under consideration at 
all relevant temperatures down to T = and at any in- 



terface transmission. 
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